TDA Progress Report 42-110 


August 15, 1992 



19 431 ,$ 


Multiple Symbol Partially Coherent Detection of MPSK 

M. K. Simon 

Telecommunications Systems Section 

D. Divsalar 

Communications Systems Research Section 


In this article , it is shown that by using the known (or estimated) value of car- 
rier tracking loop SNR in the decision metric , it is possible to improve the error 
probability performance of a partially coherent multiple phase-shift-keying (MPSK) 
system relative to that corresponding to the commonly used ideal coherent decision 
rule. Using a maximum-likelihood approach , an optimum decision metric is derived 
and shown to take the form of a weighted sum of the ideal coherent decision metric 
(i.e.j correlation) and the noncoherent decision metric which is optimum for differ- 
ential detection of MPSK. The performance of a receiver based on this optimum 
decision rule is derived and shown to provide continued improvement with increas- 
ing length of observation interval (data symbol sequence length). Unfortunately } 
increasing the observation length does not eliminate the error floor associated with 
the Unite loop SNR . Nevertheless , in the limit of infinite observation lengthy the av- 
erage error probability performance approaches the algebraic sum of the error floor 
and the performance of ideal coherent detection , i.e. } at any error probability above 
the error floor, there is no degradation due to the partial coherence. It is shown 
that this limiting behavior is virtually achievable with practical size observation 
lengths. Furthermore, the performance is quite insensitive to mismatch between 
the estimate of loop SNR (e.g. t obtained from measurement) fed to the decision 
metric and its true value. These results may be of use in low-cost Earth-orbiting or 
deep-space missions employing coded modulations. 


I. Introduction 

It is well known that for ideal phase coherent detection 
of multiple phase-shift-keying (MPSK), the decision rule 
that minimizes average bit error probability is based on 
a correlation metric and leads to bit-by-bit decisions. In 
practical situations, the phase introduced by the transmis- 
sion over the channel is unknown and thus the assumption 
of perfect knowledge of this parameter at the receiver is 
idealistic. Typically, if the channel phase is reasonably 
well behaved, the receiver will attempt to estimate it via 


some type of phase synchronization subsystem, such as a 
carrier phase tracking loop. Since the estimate is made in 
the presence of the ever-present additive channel thermal 
noise, the receiver’s phase estimate used for demodulation 
purposes is not perfect. Detection under these circum- 
stances is known as partially coherent detection. 

Ordinarily in this environment, one continues to use the 
ideal coherent detection correlation metric despite the fact 
that it is no longer optimum for partially coherent detec- 
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tion. In particular, the presence of a phase error between 
the true channel and the receiver’s estimate of it intro- 
duces memory into the observation, and thus any metric 
leading to bit-by-bit detection cannot be optimum. In- 
stead, one must resort to sequence estimation where the 
length of the sequence is proportional to the duration over 
which the phase error can be assumed constant. 


where r = (ro, v\ , nv-i), s — ( 5 o, ■■■* s Af-i)> 

n = (n 0 , n lf x ) are the received sequence, transmit- 
ted sequence, and noise sequence, respectively. Also, n k is 
a sample of zero mean complex Gaussian noise with vari- 
ance (per dimension) <7^ = No /T where No is the single- 
sided power spectral density of the noise process n(f) at 
the receiver input. 


In this article, a maximum-likelihood approach to par- 
tially coherent detection is taken, an approach not unlike 
that previously applied to noncoherent and coherent de- 
tection. It will be shown that considerable performance 
improvement can be gained by using the optimum metric 
which leads to a maximum- likelihood sequence estimation 
(MLSE) type of algorithm. 


II. Maximum-Likelihood Partially Coherent 
Detection of MPSK Over an AWGN 
Channel 

Consider the transmission of MPSK signals over an ad- 
ditive white Gaussian noise (AWGN) channel. The base- 
band representation of the transmitted signal in the inter- 
val (JfcT, (it + 1)T) has the complex form 


st = y/2 Pei*” 


( 1 ) 


where P denotes the constant signal power, T denotes 
the MPSK symbol interval, and <f> k the transmitted phase 
which takes on one of M uniformly distributed values 
j3 m = 2irm/M;m = 0, 1, . . . ,M - 1 around the unit circle. 
Assume that in addition to AWGN, the channel introduces 
a phase 0 which can be constant (independent of time) over 
a duration of N data symbols and uniformly distributed 
in the interval Thus, the received sequence r is 

expressed as 


r — se^ B + n 


( 2 ) 


For partially coherent detection, the receiver provides 
a carrier phase synchronization subsystem, e.g., a tracking 
loop, which derives a complex reference signal e* B whose 
phase 6 is an estimate of the unknown channel phase 6 . 
After demodulating r with this reference (complex conju- 
gate multiplication of the two signals), one gets 


R = Te~ j0 = se j4>c + ne j§ (3) 


where <j> c = 0 - 0 is the carrier phase error and typically 
has a Tikhonov probability density function (pdf) [1], i.e., 


P(<t>c) = 


exp(p cos <f>c) ' 

2tt/o 0) 


|4e| < » 


(4) 


Here p is a parameter related 1 to the tracking loop SNR 
and Iq(.) is the zeroth-order modified Bessel function of 
the first kind. 


For the assumed AWGN model for n, the a posteriori 
probability of the demodulated received sequence R given 
the transmitted sequence s and the carrier phase error <f> c 
follows from Eq. (3) and is 


1 For first-order tracking loops, p is indeed the loop SNR. For second- 
order loops, p is approximately the loop SNR for sufficiently large 
values [1]. In what follows, p is referred to simply as the loop SNR 
which is assumed to be known or estimated. 


p( RM«)= 1 - 2 tn ex P 

(2 *K) 


R — se J ^ c 

2^r~ 


(2ttct2) 


N 


exp < - 


2 *1 


y j + k*-»r] - 2 

.1=0 


N - 1 


y] Rk-i s k-i 


COS ( 4> c - Of) 


(5) 


i=0 
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where 



Averaging Eq. (5) over the pdf in Eq. (4) gives, upon simplification, 


( 6 ) 


p(R- I s ) = J * 


p(R\s,<t>c)p(<t>c) 





Since 


N-l 


J2 R *-i s i-i 

i—0 


cos a = Re 



( 7 ) 


( 8 ) 


Eq. (7) further simplifies to 


p(R|s) 


1 1 
h{p) ( 2 n(r%) N 



x Iq 



( 9 ) 


Note from Eq. (1) that for MPSK, |s fc | 2 is constant for all transmitted phases (3 m . Thus, since 7 0 (x) is a monotonic 
function of its argument, maximizing p(R|s) over s is equivalent to finding 
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This, using Eq. (1), results in the decision rule 


choose 


4>k-N+i if 


JV— 1 

£ Rk-ie-i*'-* 
1=0 


■ 2 ^ 


Re 



is maximum 


(ii) 


where <£k-jv+i i s a P&rficul^r sequence of the 

transmitted phases /? m . In Eq. (11), the first term in- 
side the braces represents the component of the decision 
metric associated with noncoherent (differential) detec- 
tion [2], i.e., total lack of knowledge of the uniformly dis- 
tributed channel phase 0. The second term inside the 
braces represents the component of the decision metric 
for ideal coherent detection, i.e., complete knowledge of 
the channel phase 6. Thus, the partially coherent decision 
metric is a linear combination of the coherent and non- 
coherent decision metrics with the weighting of the two 
terms in proportion to the product of the tracking loop 
SNR and the channel noise variance. Note that for any 
nonzero value of p ) this decision rule is unique because the 
second term inside the braces in Eq. (11) is unique but not 
the first term. For p = 0, which corresponds to differen- 
tially coherent detection, there is a phase ambiguity since 
the addition of an arbitrary fixed phase, say <£ a > to a h N 
estimated phases results in the same 

decision for <f>. In [2], the authors observed that by letting 
<t>a ~ ai *d differentially encoding the input phases 

at the transmitter, 

<pk = i + ( 12 ) 

where now A<£k denotes the input data phase correspond- 
ing to the ith transmission interval and 4>k the differen- 
tially encoded version of it, the decision rule can turn into 
one in which the phase ambiguity is resolved. From now 
on, assume p / 0 and thus that there is no formal require- 
ment for differentially encoding the data phase symbols. 


has perfect knowledge of p, and thus should outperform a 
conventional bit-by-bit correlation receiver which does not 
make use of this knowledge. The following sections deter- 
mine how much the optimum partially coherent sequence 
receiver outperforms the conventional bit-by-bit correla- 
tion receiver. 


III. Bit Error Probability Performance 

To obtain a simple upper bound on the average bit error 
probability, P*, of the proposed N - bit detection scheme, 
use a union bound analogous to that used for upper bound- 
ing the performance of error correction coded systems. 
In particular, the upper bound on is the sum of the 
pairwise error probabilities associated with each AT-bit er- 
ror sequence. Each pairwise error probability is then ei- 
ther evaluated directly or itself upper bounded. Math- 
ematically speaking, let = (</>k,<f>k-u • • • »^*-~JV+i) de- 
note the sequence of N transmitted information phases 
and (^,^-i,...,^*-iv+i) be the corresponding se- 
quence of detected phases. Let u be the sequence of 
6 - N log 2 M information bits that produces <f> at the 
transmitter and let u be the sequence of b bits that results 
from the detection of $ Then, since MPSK is a symmetric 
signalling set, i.e., it satisfies a uniform error probability 
(UEP) criterion, one gets an upper bound on the bit error 
probability, 

< N lo g 2 ^ 7 J2 w ( u - & ) Pr {v>n\4>,<t>c} (13) 


Figure 1 is an illustration in complex form of a re- 
ceiver implemented on the basis of Eq. (11). Note that 
this receiver requires knowledge of the loop SNR p } the 
signal power P, and the noise variance cr 2 . The accuracy 
of this knowledge, which must be obtained by measure- 
ment, will have an impact on the ultimate performance 
of this receiver. Later, in Subsection E, the authors in- 
vestigate the sensitivity of the receiver to a mismatch be- 
tween the true loop SNR and the value supplied to the 
receiver implementation in Fig. 1. In the next section, 
except in Section III.E, it is assumed that the receiver 


where the decision statistic tj is defined from Eqs. (10) and 
(11) by 2 


V = 


JV-1 






1 = 0 



(14) 


2 Note that when compared with Eq. (11), v of Eq. (14) include* the 
additional constant (pcr^/v^P) 2 . This, however, has no effect on 
the decision-making process and thus one can use the convenient 
form of Eq. (14) with no loss in generality. 
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and the corresponding error statistic fj is identical to 
Eq. (14) with each 0* replaced by In Eq. (13), w(u,u) 
denotes the Hamming distance between u and u, <f> is 
any input sequence (e.g., the null sequence (0,0,..., 0) = 
0), and Pr {j 7 > r; | denotes the pairwise probabil- 


ity that <f> is incorrectly chosen when indeed 4> was 
sent. Note that the bound in Eq. (13) is computed for 
a fixed carrier phase error, <f> Ci which accounts for the 
notational dependence of Pr {77 > 77 | <j> y <t> c ) and thus 
P b (<j>c) on <£ c . 


A. Evaluation of the Pairwise Error Probability 

To compute Pr {rj > 77 | <£, <£ c }, the approach taken in [2] is used for evaluating the performance of multiple symbol 
differentially coherent detection of MPSK. In particular, letting i] = |zi| 2 and f) = |*2| 2 [see Eq. (14) and the statement 
below it for the definitions of z\ and z 2], then [3] 

Pr {17 > r\ | <£,<£«.} = - |l - Q{Vb ) >fa) + Q(\/a, \/t)j = f(a,b) (15) 


where Q(x,y) is the Marcum Q function [4] and 


f 6 | _ 1 (5x^52-21^17? ^ cos (^-ga-ft/) 5i - 5 2 | 


where the + sign and — sign correspond to b and a, respectively, and 


Si 


N + 


2E,/N 0 


0-11 


= p 



p 

E,/N 0 


N COS 0 e 


+ 



(16) 


S 2 = P 


5 + 


2 E,/N 0 




= P |fi| 8 + 


E./N 0 


Re 


{6^} + ( 


2E,/N 0i 


li =2 m N 0 

Nz = 2 l Zl “ 2 il = N ~fr 

t - rfirWi - ( 2 2 - W = ^ = ar 8 ( — arg 6 

9 ,=argTr=»g (Ar e «. + _£_| ; S, = „g = ..g (17) 


and 


N—l 

8 = ^ ei(^-.-^-Oo (18) 

i=0 

which is a normalized time cross correlation between s and s. Also, E,/No — PT/No denotes the symbol energy-to-noise 
spectral density ratio and is related to the bit energy-to-noise spectral density ratio Eh/No by E,/Nq = (Eb/No) log 2 M. 
Substituting Eq. (17) into Eq. (16) results, after considerable simplification, in 
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!PP[ 


{*.}' -k H‘- 4 w- -* + W<'-H} 

I El fZ 3 | f | a I I 1 ( - — ) (jtfcoa^e - |g|cos(^ + y)) ( 19 ) 

± m[v K w \jip-\tf w" J 

Now some special cases of practical interest are considered. 

B. Case 1: Binary PSK With Two-Symbol Observation and Detection (M= 2, W = 2) 

In this case, E./N 0 = E b /N 0 . There are M 2 - 1 = 3 possible error sequences each of length 2. The pertinent results 


related to the evaluation of Eqs. (18) and (19) are 


<j>k - 4>k 4>k-\ - 4>k-i & 


For the first two error sequences, Eq. (19) evaluates to 


b =£k 4+2 


(m) w ‘ + K 


p 

E b /N 0 


4 \E b /NoJ 


For the third error sequence, both a and 6 approach infinity (the ratio a/6, however, approaches unity) as 6 
-2. Thus, one must evaluate the pairwise error probability Eq. (15) separately for this case. I is g 

show that 


lim f(a, 6) — ^ erfc 


a — 00 
b—+OG 
a/6-1 


Furthermore, in the general case where 6 -* -N , Eq. (21) evaluates to 




which for N = 2 and M = 2 becomes 


lim /(a, 6) = - erfc 
6—* — 2 « 
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Finally, noting that the Hamming distance w(u,u) is equal to 1 for the first two error sequences and is equal to 2 for 
the third sequence, substituting Eqs. (23) and (15) combined with Eqs. (19) and (20) into the expression for bit error 
probability in Eq. (13) gives 


PM < f 


(- 

\N o 


8 ( £»/*«,) 


E b _ 
'No 




COS <j> c + - 



+ - erfc 



(24) 


Finally, the upper bound on average bit error probability Pi, is obtained by averaging the upper bound in Eq. (24) over 
the pdf in Eq. (4). Figures 2 and 3 are plots of this upper bound on average probability versus E^/Nq in decibels for 
values of p — 7 dB and 10 dB, respectively. For the purpose of comparison, the exact results (i.e., not an upper bound) 
for the conventional ideal coherent metric operating in a noisy carrier synchronization environment are 


Pi 




p(<j> e )d<j>c 


(25) 


where p(<f> c ) is given by Eq. (4). Even with only one additional observation symbol interval, considerable savings in 
Eb/No can be achieved at a fixed error probability, particularly in the region of the knee of the curve where the system 
begins approaching its irreducible error probability 3 asymptote (error floor). 

C. Case 2: Binary PSK with W-Symbol Observation and Detection (M= 2 , N arbitrary) 

For N arbitrary, 6 takes on values -(N - 2 i)\i - 0,1, 2, ...,1V - 1. The number of error sequences corresponding 
to each of these values of 6 is binomially distributed, i.e., there are (^) sequences that yield a value 6 = —(N — 2«). 

Furthermore, the Hamming weight associated with each of the (^) sequences that yield a value 6 = —(N — 2 *) is 
tv(u,u) = N — i. Finally then, using the above in Eqs. (19) and (22) and substituting the results in Eq. (13), the 
conditional bit error probability is upper bounded by 


We) < 


^ erfc (\f^¥ 0 cos4> ) + e (f)(^- i )/(«.-.i.) 

= ^ erfc ^ ^ ^ /(a^t.) 


(26) 


where 


{ £ } - m, { [ N + tak) c “ ' *• + i (mx) ] * + ( im) ' * •] } (27) 


3 It is well known [1] that conventional PSK systems exhibit an irreducible error probability (i.e., a finite error probability in the limit 
as E^/No approaches infinity) when a noisy carrier synchronization reference with fixed power is used as a demodulation signal. This is 
observed by examining a curve of P versus E^/No with loop SNR, p, held fixed. The value of this irreducible error probability is given 
by [1] Pblirr = /// 3 p(^c)<tyc Note that in practice, as the observation length increases, one should decrease the loop bandwidth of the 
phase-locked loop (PLL), which results in an increase in the loop SNR, Also, as the bit SNR increases, the loop SNR (for fixed modulation 
index) increases and thus the error floor decreases. 
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The upper bound on unconditional average bit error probability is now obtained by averaging Eq. (26) over the pdf in 
Eq. (4). The numerical results are illustrated in Figs. 2 and 3 for values of N = 4, 6, and 8. As N gets large, the curves 
appear to approach an asymptote. This asymptotic behavior is analytically evaluated as follows: 

For large N, the first term in Eq. (26) when integrated over the pdf in Eq. (4) approaches the irreducible error 
probability P b \ irr as f* /a p(4> c )d<j> e . Also, the dominant term in the summation term of Eq. (26) corresponds to i = N-l, 
i.e., S = N - 2. Thus, for large N, the second term of Eq. (26) approaches /(ajv-i,&w-i) where 

{t»:‘,} e! ^{ w * 2 ' /wri >} (28) 

Since from Eq. (28), y/hfl 7 » y/b^-i - VajvTT, then using the asymptotic form of Eq. (15) for a and b large (see 
Appendix A of [5]), namely, 

/(«.»)* 5 < 29 > 

The value 5 t^-i) w obtained as 

f(a N -i,b N .i) 2 ^ erfc ( 30 ) 

independent of <f> e . Finally then, for large N, the asymptotic behavior of the average bit error probability is approximately 
upper bounded by 

(31) 

namely, the mm of the bit error probability for ideal coherent detection and the error floor. Equation (31) is in very close 
agreement with the curves for N = 8 in Figs. 2 and 3. 

D. Case 3: Quaternary PSK With TVvo-Symbol Observation and Detection {M = 4, N = 2) 

In this case, E,/N 0 = 2E b /N 0 - There are now a total of M 2 - 1 = 15 possible error sequences each of length 2. Of 
these, only eight produce distinct combinations of |6| and v. These are tabulated below: 


Error 

sequence 

Case 

<t>k ~ 

d>k-i — 4>t-i 

\6\ 

V 

1,2, 3, 4 

i 

;r,0, 3 jt/ 2, ir/2 

0, ir, jt/2, 3it/2 

0 

0 

5,6 

2 

0,ir/ 2 

n/2,0 

\/2 

?r/4 

7,8 

3 

0,3tt/2 

3tt/2,0 

n/2 

-tt/4 

9,10 

4 


5r, jt/2 

V2 

3*r/4 

11,12 

5 

7 r, 3 jt/2 

3tt/2, tt 

V2 

— 3tt/4 

13 

6 

7 c / 2 

tt/2 

2 

7t/2 

14 

7 

3tt/2 

3ir/2 

2 

— jt/ 2 

15 

8 

7 T 

7 T 

2 

w 
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The corresponding values of a and b for each of the first five cases which correspond to |5| ^ 2 are given as follows. 


m 

Wl _ *0 




a 


{i}.- 


3 ^0 


2 +( 

P \ 

COS 4>c + i I 

p 

)'M 

^ 2 E b /N 0 ) 

<2E b /N , 0 

2 + 

( P \ 

cos <j> c + ^ 1 

f P 

VI 

4- 

K^Ei/No) 

K 2E b jN 0 

>\ 


2 + 

( P \ 

. 1 
COS <j> c + - 

f p 

yi 

± [ 

\2 E b /No) 

\2E b /N 0 

)\ 

± L 

2 + 1 

f P ) 

3 , 

| COS 4> c + - 1 

( P 

v 

± [ 

\2E b /N 0 J 

\2E b /N 0 

] . 

L 

2+1 

f P ] 

3 . 

| COS 4> c + - | 

( P 

v 

4 . 

\2E b /N 0 J 

\2E b /No 

J . 



,2 E b /N 0 


) COS 4>c | 


)1 * h + ^(^) (co * , ‘“" i "*' ) ]} 


JEi/No 


)(3 


cos <f> c - sin <f> c ) 


(32) 


For case 9 6 and 7, the following is analogous to Eq. (22); 


. l>m v /( a . & ) = o erfc 

6-*±jN 2. 


I NE„ 

2N 0 


(cos <t> c ± sin <j > c ) 


) 


(33) 


which for N = 2 and A/ = 4 becomes 


lim /(a,6) = - erfc 
<5->±j2 Z 


(y^- (cos 4> c ± sin tf c )^ 


(34) 


Finally, for case 8, Eq. (22) is used to obtain 


lim /(a, 6) = 

o* - * 2 


- erfc 
2 



(35) 


Evaluating the Hamming distances for the 15 error sequences and substituting the above results into the expression for 
the bit error probability bound in Eq. (13) gives 


P b (<t> e ) (6/ (a Xl 50 + 2/ (a 2 , 6 2 ) + 2/(a 3l fc 3 ) + 4/ (a 4) t 4 ) + 4/(a 5 A)} 



Ftote* 4 and 5 are comparable to Figs. 2 and 3 for the M = 4 (QPSK) case. The analytical exact result corresponding 
to the ideal coherent metric operating in a noisy carrier synchronization environment is now [1] 

Pt= f ierfcfy^(cos* c + »in*,)^W.W,4. £J« rfc (^(cos*« - si„<w) <*0 

Analysis and plots for larger values of N are not included here, but they would show further improvement as was true 
for the binary case. 


E. Performance Sensitivity to Mismatch 
In Loop SNR 

Here the authors investigate the sensitivity of the aver- 
age bit error probability (in terms of its upper bound) of 
the MLSE receiver to a mismatch between the true loop 
SNR, p, and the estimate of it, p, supplied to the imple- 
mentation of Fig. 1. In particular, the authors evaluate, 
for the special cases of Sections III -B and III.C, the upper 
bound 

Pb< j Pbu{4>c\p)p(<t>c)d<i>c ( 38 ) 


where Pb u {4>c,p) is given by the upper bound in Eq. (24) 

or Eq. (26) with p replaced by p = p[l + (p - p)/p] = 
p(l + e) and p(<t> e ) is as given by Eq. (4). Figures 6 and 
7 are illustrations of Eq. (38) for M = 2,p = 10 dB, and 
N — 2 and 8, respectively, with fractional mismatch e as 
a parameter. One observes that even with mismatches as 
much as 50 percent (e = ±0.5), there is negligible effect 
on the error probability performance. Thus, the authors 


conclude that the MLSE receiver is quite insensitive to 
mismatch in the loop SNR. 


IV. Conclusions 

By making use of the known (or estimated) value of 
loop SNR in the decision metric, it is possible to im- 
prove the error probability performance of a partially co- 
herent MPSK system relative to that corresponding to 
the commonly used ideal coherent decision rule. Using 
a maximum-likelihood approach, an optimum decision 
metric was derived and shown to take the form of a 
weighted sum of the ideal coherent decision metric (i.e., 
correlation) and the noncoherent decision metric previ- 
ously shown to be optimum for differential detection of 
MPSK. The performance of a receiver based on this op- 
timum decision rule improves with the increasing length 
of the observation interval (data symbol sequence length). 
Furthermore, the performance is quite insensitive to mis- 
match between the estimate of loop SNR (e.g., obtained 
from measurement) fed to the decision metric relative and 
its true value. 
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UPPER BOUND ON AVERAGE BIT ERROR PROBABILITY 







UPPER BOUND ON AVERAGE BIT ERROR P 



Fig. 4. Upper bound on average bit error probability versus E b IN 0 
In decibels for MLSE and comparison with exact performance of 
Ideal coherent metric; M= 4, Ns 2, and p - 13 dB. 
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Fig. 5. Upper bound on average bit error probability versus E b IN 0 
In decibels for MLSE and comparison with exact performance of 
ideal coherent metric; M « 4, Ns 2, and p e 16 dB. 
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UPPER BOUND ON AVERAGE BIT ERROR PROBABILITY 



5 6 7 8 9 10 11 12 13 14 15 


E b iN 0 ,dB 

Fig. 6. Upper bound on average bit error probability veraua E b IN 0 
In decibels for MUSE in the presence of loop SNR mismatch; M ■ 2, 
Ws 2, and p* 10. 



Fig. 7. Upper bound on average bit error probability versus E b /N 0 
In decibels for MLSE In the presence of loop SNR mismatch; M s 2, 
N = 8, andpslO. 
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